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ABSTRACT 


This  report  describes  ground  moving  target  indicator  (GMTI)  processing  to  be  done  for  the  PCA 
(polymorphous  computing  architecture)  integrated  radar-tracker  application.  GMTI  processing  of  the  raw 
radar  data  is  done  to  extract  targets;  target  reports  are  passed  on  to  the  tracker  part  of  the  integrated  radar- 
tracker  application. 

The  GMTI  processing  described  in  this  report  is  used  to  minimize  the  effects  of  interference  patterns 
and  ground  clutter  and  to  recognize  and  extract  targets  on  the  ground.  The  radar  is  presumed  to  be  using  a 
wideband  signal  for  improved  resolution.  This  signal  is  decomposed  via  subband  filtering  into  many 
narrowband  signals  for  easier  processing.  Each  of  these  subbanded  data  sets  must  then  be  processed  as  a 
distinct  radar  data  cube.  For  each  data  set,  adaptive  processing  is  used  to  null  interference  patterns  and 
ground  clutter  noise.  These  operations  incur  a  high  computational  cost.  Subband  filtering  requires  many 
parallel  modulators  and  filters,  while  adaptive  processing  requires  the  calculation  of  an  adaptive  filter 
(involving  a  QR  decomposition  and  forward/backsolving)  at  run  time  for  each  subband  data  cube. 
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1.  INTRODUCTION 


This  report  describes  ground  moving  target  indicator  (GMTI)  processing  to  be  done  for  the  PCA 
(polymorphous  computing  architecture)  integrated  radar-tracker  application.  GMTI  processing  of  the  raw 
radar  data  is  done  to  extract  targets;  target  reports  are  passed  on  to  the  tracker  part  of  the  integrated  radar- 
tracker  application. 

The  GMTI  processing  described  in  this  report  is  used  to  minimize  the  effects  of  interference  patterns 
and  ground  clutter  and  to  recognize  and  extract  targets  on  the  ground.  The  radar  is  presumed  to  be  using  a 
wideband  signal  for  improved  resolution.  This  signal  is  decomposed  via  subband  filtering  into  many 
narrowband  signals  for  easier  processing.^  Each  of  these  subbanded  data  sets  must  then  be  processed  as  a 
distinct  radar  data  cube.  For  each  data  set,  adaptive  processing  is  used  to  null  interference  patterns  and 
ground  clutter  noise.  These  operations  incur  a  high  computational  cost.  Subband  filtering  requires  many 
parallel  modulators  and  filters,  while  adaptive  processing  requires  the  calculation  of  an  adaptive  filter 
(involving  a  QR  decomposition  and  forward/backsolving)  at  run  time  for  each  subband  data  cube. 

This  report  describes  the  key  features  and  elements  for  a  preliminary  design  review  of  the  GMTI 
radar  processing  component  of  the  integrated  moving  target  indicator/tracker  application.  This  report 
includes: 

•  a  high-level  description  of  the  GMTI  radar  processing  chain 

•  specifications  of  the  input  and  output  for  each  of  the  GMTI  radar  processing  chain  elements 

•  detailed  descriptions  of  each  of  the  GMTI  radar  processing  chain  elements 

1.1  THE  GMTI  RADAR  PROCESSING  CHAIN 

Before  delving  into  the  details  of  each  of  the  elements  of  the  processing  chain,  a  description  of  the 
entire  GMTI  processing  chain  is  in  order.  The  processing  chain  is  illustrated  in  Figure  1 .  In  this  figure,  the 
small  numbers  in  the  upper  right  comer  of  the  boxes  refer  to  the  section  in  which  that  stage  is  discussed  in 
greater  detail.  The  shaded  box  with  bold  labels  is  described  in  [1]. 


*  A  wideband  signal  is  a  signal  in  which  the  bandwidth  is  a  significant  fraction  of  the  center  frequency,  while  a 
narrowband  signal  is  one  in  which  the  bandwidth  is  an  insignificant  fraction  of  the  center  frequency.  The  determination 
of  what  constitutes  a  significant  fraction  of  the  center  frequency  is  application  specific.  As  this  is  a  wideband  radar, 
the  assumptions  made  for  narrowband  processing  break  down,  leading  to  dispersion  in  the  beams  that  are  formed. 
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Figure  I.  GMTI processing  chain. 


Before  the  radar  data  is  processed  by  the  GMTI  stages,  a  radar  signal  consisting  of  a  series  of  pulses 
comprising  a  coherent  processing  interval  (CPI)  has  been  transmitted.  The  pulse  repetition  interval  (PRI) 
determines  the  time  interval  between  transmitted  pulses,  and  the  number  of  pulses  per  CPI  is  given  as  . 
The  signal  reflects  off  of  targets,  the  ground,  and  other  entities  and  is  received  by  the  radar  antenna.  The 
antenna  is  an  array  of  sensors,  which  are  also  called  channels.  These  channel  sensor  inputs  are 
digitized  by  A/D  converters  at  a  rate  of  samples  per  PRI.  The  digitization  rate  alludes  to  the 
relationships  between  the  range  dimension  and  the  PRI  dimension:  the  range  dimension  is  referred  to  as  the 
fast  time  dimension,  while  the  PRI  dimension  is  the  slow  time  dimension.  This  digitized 

data  cube  of  samples  is  the  input  to  the  GMTI  processing  chain. 

The  GMTI  processing  chain  consists  of  nine  stages:  subband  analysis  (decomposing);  time  delay 
and  equalization;  adaptive  beamforming;  pulse  compression;  Doppler  filtering;  space-time  adaptive 
processing  (STAP);  subband  synthesis  (recombining);  detection;  and  estimation.  This  processing  chain 
then  reports  the  resulting  targets  to  the  tracker. 

•  The  suhband  analysis  stage  breaks  the  wideband  radar  signal  into  a  series  of  narrowband  signals 
to  prevent  beam  dispersion  which  facilitates  more  efficient  processing. 

•  The  time  delay  and  equalization  stage  compensates  for  differences  in  the  transfer  function 
between  channel  sensors. 
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•  The  adaptive  beamforming  stage  transforms  the  filtered  data  into  the  beam-space  domain  to 
allow  detection  of  target  signals  coming  from  a  particular  set  of  directions  of  interest  while 
filtering  out  spatially-localized  interference. 

•  The  pulse  compression  stage  filters  the  data  to  concentrate  the  signal  energy  of  a  relatively  long 
transmitted  radar  pulse  into  a  relatively  short  pulse  response. 

•  The  Doppler  filter  stage  processes  the  data  so  that  the  radial  velocity  of  targets  relative  to  the 
platform  can  be  determined. 

•  The  STAP  stage  is  a  second  beamforming  stage  which  removes  further  interference  and  ground 
clutter  interference. 

•  The  subband  synthesis  stage  stitches  the  processed  subband  signals  back  together,  reversing  the 
work  that  the  subband  analysis  stage  had  conducted  on  the  unprocessed  signal. 

•  The  detection  stage  uses  constant  false-alarm  rate  (CFAR)  detection  to  compare  a  radar  signal 
response  to  its  surrounding  signal  responses  to  determine  whether  a  target  is  present  and  uses 
target  grouping  to  eliminate  multiple  target  reports  that  are  actually  just  one  target. 

•  The  estimation  stage  estimates  target  positions  in  order  to  pass  them  to  the  tracking  algorithms. 

•  Finally,  the  target  report  is  passed  to  the  tracking  algorithms  which  are  described  in  [1]. 

1.2  GMTI  CONTROL  INPUT  DATA 

The  GMTI  processing  function  requires  a  control  input  which  is  a  Matlab  structure.  The  names  of 
the  elements,  the  variables  they  correspond  to  in  this  report,  and  the  descriptions  of  each  element  are  given 
in  Table  1 .  This  structure  must  be  passed  to  the  GMTI  fiinction  when  it  is  called. 
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TABLE  1 


Control  Input  Parameters  for  GMTI  Processing 


Parameter 

Element 

Name 

Description 

Nch 

Number  of  radar  sensor  channels 

Nrg 

Number  of  range  gates  per  radar  data  cube 

Npri 

Number  of  pulse  repetition  intervals  per  radar  data 
cube 

Nsub 

Number  of  subbands  for  processing 

Nsrg 

Number  of  subband  range  gates  in  radar  data  cube 

Ndwn 

Down-sample  rate  for  the  subband  analysis  stage,  and 
conversely,  the  up-sample  rate  for  the  subband 
synthesis  stage 

^sfd 

Nsfd 

Number  of  finite  impulse  response  (FIR)  filter  taps  in 
the  subband  analysis  down-sampling  low-pass  filter 

Nsfu 

Nsfu 

Number  of  FIR  filter  taps  in  the  subband  synthesis  up- 
sampling  low-pass  filter 

N,de 

Ntde 

Number  of  FIR  filter  taps  in  the  time  delay  and 
equalization  filter 

a 

alpha 

Diagonal  loading  factor  in  the  adaptive  beamforming 
stage 

Npc 

Number  of  filter  taps  in  the  pulse  compression  FIR 
filter 

f^bn, 

Nbm 

Number  of  beams  in  the  radar  data  cube  after  the 
adaptive  beamforming  stage 

Ndop 

Ndop 

Number  of  Doppler  bins  for  Doppler  filtering 

^stag 

Number  of  PRI  staggers  produced  by  Doppler  filtering 
for  use  in  the  STAR  stage 

N,s 

Number  of  training  samples  for  the  STAR  stage 

Knb 

Ncnb 

Number  of  clutter  nulled  beams  in  the  radar  data  cube 
after  the  STAR  stage 

^cfar 

Ncfar 

Number  of  range  gates  used  in  forming  the  noise 
estimate  from  each  side  of  the  cell  under  test 

G 

G 

Number  of  guard  range  cells  on  both  sides  of  the  cell 
under  test 

M 

mu 

Normalized  target  power  threshold 

Most  of  the  control  input  parameters  presented  in  the  table  above  are  also  shown  in  tables  at  the  end 
of  the  stage  section  in  which  they  are  used. 

Beyond  the  control  input  parameters,  there  are  a  number  of  data  product  inputs,  as  shown  in  Table  2. 


TABLE  2 

Data  Product  Inputs  for  GMTI  Processing 


Parameter 

Element 

Name 

Description 

hs(n) 

hs 

Coefficient  vector  of  subband  analysis  low-pass  filter 

fs^n) 

f  s 

Coefficient  vector  of  subband  synthesis  low-pass  filter 

'^bm 

Tbm 

Vector  of  range  gate  indices  for  adaptive  beamforming 
stage  training  samples 

\c 

hpc 

Coefficient  vector  of  the  pulse  compression  filter 

'^stag 

Tstag 

Vector  of  space-time  snapshot  indices  for  the  training 
data  for  space-time  adaptive  processing 

13  GMTI  INPUT  DATA 

The  input  data  to  the  GMTI  radar  processing  chain  is  a  series  of  three-dimensional  radar  data  cubes 
whose  dimensions  are  channel,  range,  and  pulses  (see  Figure  2).  The  collection  of  the  samples  of  the  radar 
data  cubes  is  explained  in  Section  1.1.  Individual  elements  of  the  radar  data  cube  C  are  referred  to  as  C(i,j, 
k),  where  i  is  the  channel  index,  y  is  the  range  index,  and  k  is  the  pulse  index. 

The  series  of  radar  data  cubes  will  be  indexed  using  a  cell  array  in  Matlab: 
input_data { 1 }  =  cubel ; 
input_data{2 }  =  cube2 ; 
input_data{ 3 }  =  cube3 ; 
etc . 


^  Items  in  Courier  font  correspond  to  functions  or  variables  in  the  Matlab  code. 
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Figure  2.  Radar  data  cube,  the  input  data  to  the  GMTI  processing. 


1.4  GMTI  OUTPUT  DATA 

The  tracker  receives  the  following  data  on  each  target: 

•  target  power  or  signal-to-noise  ratio  (SNR) 

•  target  azimuth  (in  degrees  or  radians) 

•  target  range  (in  meters  or  kilometers) 

•  target  radial  velocity  (in  m/s  or  km/h) 

•  target  time  stamp 

The  tracker  is  described  in  [1], 
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2.  FUNCTIONAL  DESCRIPTION 


In  this  section,  each  of  the  processing  stages  are  described  in  greater  detail.  For  convenience,  Figure 
3  shows  the  GMTI  processing  chain  with  the  sizes  of  the  passed  parameters. 


[g 
— 1^ 

Compute 

STAR 

- 

Weights 

Figure  3.  GMTI  processing  chain  with  passed  parameter  sizes. 


2.1  SUBBAND  ANALYSIS  (AND  SYNTHESIS) 

Input:  Wideband  radar  data  cube  of  size  x  . 

Output:  Narrowband  subband  data  cubes  of  size  x 

In  order  to  get  greater  resolution  in  recognizing  targets,  it  is  beneficial  to  use  wideband  radar  pulses. 
However,  channel  equalization,  pulse  compression,  Doppler  filtering,  and  adaptive  nulling/beamforming 
must  be  performed  on  a  narrowband  signal  because  the  assumptions  made  for  narrowband  processing 
break  down  in  a  wideband  scenario,  leading  to  dispersion  in  the  beams  that  are  formed.  (Compensation  for 
beam  dispersion  makes  a  system  far  more  complex.)  Therefore,  the  processing  chain  must  involve 
decomposing  the  wideband  signal  into  many  downsampled  narrowband  signals,  processing  each  subband 
individually,  and  then  recombining  the  subband  signals  prior  to  performing  target  detection  [6].  The 
subband  filter  is  implemented  with  a  classic  polyphase  filter  [5].  Since  subband  analysis  and  synthesis  are 
duals  of  each  other,  both  stages  will  be  discussed  within  this  section. 
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Subband  analysis  consists  of  modulating  to  baseband  the  incoming  signal  into  a  series  of  subband 
signals,  low-pass  filtering  each  of  these  subband  signals  with  h^(n) — an  Af^y^-tap  filter — and 
downsampling  the  subsequent  filtered  signals.  The  inputs  to  this  function  are  the  range  vectors  of  the  input 
radar  data  cube;  x  =  C(i,  k),  where  x(n)  =  C(i,  n,  k).  The  modulation  is  performed  by  multiplying  each 
sample  by  a  fi-equency  down-shifting  value  where  and  The  low-pass 

filtering  is  conducted  so  that  signal  aliasing  is  minimized  in  the  down-sampling  step.  Down-sampling  is 
conducted  by  simply  extracting  one  out  of  every  „  samples.  When  these  operations  are  completed,  the 
output  is  a  set  of  subband  data  cube  of  size  where  N^rg  =  Kg/Nj^.„.  This 

transformation  is  depicted  in  Figure  4.  In  Matlab  notation,  these  data  cubes  will  be  referred  to  as 
k,  s),  where  /  is  the  channel  index,  y  is  the  subband  range  index,  k  is  the  PRI  index,  and  s  is  the 
subband  index.  Often  and  are  chosen  such  that  this  means  that  the  total  number  of 

samples  is  preserved  between  the  input  data  cube  and  the  subband  data  cubes. 


Subband  1  Subband  Ngub 


Figure  4.  Data  size  transformation  in  suhband  analysis  stage. 


Conversely,  subband  synthesis  involves  the  dual:  up-sampling  the  incoming  subband  signals,  low- 
pass  filtering  the  up-sampled  subband  signals  with  //«)— an  N-tap  filter— modulating  the  subband 

signals  back  up  to  their  original  frequencies,  and  recombining  the  subbands  by  superposition  addition. 
Upsampling  is  performed  by  inserting  - 1  zero-valued  samples  between  each  sample.  Again,  low-pass 
filtering  is  conducted  to  minimize  signal  aliasing  in  the  up-sampling  step.  The  modulation  is  performed  by 
multiplying  each  sample  by  a  frequency  up-shifting  value  ,  where  Wfj  =  and  0  <  it  < 

The  subband  analysis  and  synthesis  processes  are  depicted  in  Figure  5. 
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Figure  5.  GMTI  subband  filtering  stage  process. 


This  structure  is  required  for  each  channel  of  the  system.  There  is  a  great  deal  of  redundancy  within 
each  channel’s  subband  processing,  especially  in  the  low-pass  filters.  Each  of  the  h/n)  downsampling  and 
f/n)  upsampling  low-pass  filters  have  the  same  coefficients.  Therefore,  some  novel  signal  processing 
“tricks”  can  be  implemented  to  reduce  the  operations  count  [7],[8].  However,  these  methods  are  beyond 
the  scope  of  this  report. 

The  data  parameters  of  the  subband  filtering  stages  are  listed  in  Table  3. 


TABLE  3 

Data  Parameters  of  the  Subband  Filtering  Stages 


Name 

Description 

Nsut 

Number  of  subbands  for  processing 

Nsfd 

Number  of  taps  for  subband  down-sampling  low-pass  filter 

Nsfu 

Number  of  taps  for  subband  up-sampling  low-pass  filter 
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2.2  TIME  DELAY  AND  EQUALIZATION 


Input:  subband  data  cubes  of  size  x  x  . 

Output:  subband  data  cubes  of  size  x  x  . 

The  time  delay  and  equalization  stage  is  applied  to  each  subband,  and  it  uses  a  finite  impulse 
response  (FIR)  filter  on  each  range  vector  for  each  channel,  each  pulse,  and  each  subband  to  compensate 
for  signal  diflFerences  between  channel  sensors.  That  is,  in  Matlab  notation,  it  operates  on  each  x  =  C(i, k, 
s),  where  i  is  the  channel  index,  k  is  the  PRI  index,  and  5  is  the  subband  index.  In  its  simplest  form,  a  single 
tap  adjustment  can  be  folded  into  the  digital  beamforming  weights;  however,  more  complex  equalization 
requires  using  programmable  FIR  filters.  The  FIR  filter  has  filter  coefficients,  *,^4/),  1.  The 

output  of  the  filter,  the  vector  j;,  is  the  convolution  of  the  filter  with  the  input  x : 

y\J]=  =  0, 1, 

/  =  o 


Depending  on  the  number  of  filter  coefficients,  this  filter  can  be  implemented  using  fast 
convolution  with  FFTs  [4]  or  by  explicitly  convolving  the  input  signal  and  the  filter  coefficients  as 
depicted  in  the  equation  above. 

The  time  delay  and  equalization  stage  parameter  is  listed  in  Table  4. 

TABLE  4 


Time  Delay  and  Equalization  Stage  Parameter 


Name 

Description 

N,de 

Number  of  taps  in  the  time  delay  and  equalization  filters 

13  ADAPTIVE  BEAMFORMING 

Input:  subband  data  cubes  of  size  x  x  . 

Output:  subband  data  cubes  of  size  x  x  N^^.. 

The  adaptive  beamforming  stage  transforms  the  filtered  data  into  the  beam-space  domain  to  allow 
detection  of  target  signals  coming  from  a  particular  set  of  directions  of  interest  while  filtering  out  jamming 
(spatially-localized)  interference.  The  filtering  is  executed  for  each  subband  radar  data  cube  within  each 
CPI. 


Again,  we  follow  through  the  execution  of  this  stage  for  just  one  subband  radar  data  cube,  namely 
the  5th  subband.  To  perform  adaptive  beamforming,  a  clutter-free  data  matrix  is  extracted  from 
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Typically,  the  first  PRI  signal  of  a  CPI  is  a  receive-only  PRI;  that  is,  no  signal  is  transmitted  for  the  first 
PRI,  and  only  the  jamming  signals  are  sensed  by  the  radar  receiver.  This  means  using  the  channel  columns 
samples  from  the  range  and  pulse  set;  in  Matlab  terms,  it  is  C(:,j,  k,  s)  for  a  given  set  of  j’s  and  k=l.A  total 
of  channel  column  samples  are  taken,  and  the  choice  of  which  range  samples  to  use  is  given  by  the 
-length  index  vector  .  Thus,  is  an  x  matrix. 

In  adaptive  processing,  the  goal  is  to  compute  columns  of  the  beamforming  matrix,  (of  size 
Nc^xNi,„),  by  solving  where  is  the  w-th  column  of  is  the  estimated 

covariance  matrix,  and  is  the  ?K-th  column  of  the  steering  matrix,  To  this  end,  the  estimated 
covariance  matrix,  r!''\  is  computed  from  the  clutter-free  data  matrix  added  to  the  diagonal  loading  factor, 
a,  multiplied  by  an  x  identity  matrix.  is  then  a  x  matrix  as  shown  in  (2). 


5N, 


ch 


The  /M-th  column  of  the  steering  matrix,  is  computed  as 


(2) 


=  [*0 


e  ,  bye 


(A/c/,-ir 


(3) 


where  ^  and  ^  are  provided  for  each  subband  and  each  and  the  by  weights  are  static  and 
determined  as  part  of  the  radar  design.  The  range  for  m  is  0  <  w  <  so  is  an  Nf,^  x  matrix. 

Now  the  w-th  column  of  the  beamforming,  interference-nulling  matrix  is  computed  as: 


(•s) 

(R  )  % 

D(s)r' 


(4) 


The  denominator  of  Equation  (2)  is  a  scalar  normalization  factor.  The  above  equation  is  referred  to 
as  a  Wiener-Hopf  equation  for  calculating  adaptive  filter  weights.  While  this  explains  the  mathematics 
behind  calculating  the  interference-nulling  matrix,  calculating  the  inverse  of  .R  ^  is  not  used  due  to 
numerical  stability  issues.  Instead  LQ-decomposition  and  forward-  and  back-substitution  is  used  to 
calculate  each  of  the  vectors. 

As  explained  in  Appendix  A,  define  x' =  [x^/7a]>  ^  ^  by  lx67V^;,  vector.  With  (x')(y)^  =  R, 
V  =  lqhouse{{,x'){x'f)',  V  is  a  lower  triangular  N^,,xl^ch  matrix.  Now  using  as  a  temporary  vector  of 
size  N^i,  by  1,  we  can  solve  L't^  =  for  using  forward  substitution  and  for  using 

back-substitution.  Finally,  compute  take  into  account  the  denominator  scalar 

normalization  factor. 

Refer  to  Appendix  A  for  a  deeper  explanation  of  this  method  of  solving  Wiener-Hopf  equation. 
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Finally,  is  used  to  beamform  the  radar  data  cube.  Using  Matlab  notation  in  defining  the  matrix 
4*'  =  C(:, k,  s)  with  size  the  matrix  as  the  output  of  the  beamforming, 


(5) 


Thus,  the  matrix  is  of  size  x  .  Figure  6  shows  a  conceptual  data  cube  and  its  partitioning 
for  this  operation. 


N^„  Channels 


Subband  s 

Figure  6.  Data  partitioning  of  input  data  cube  for  adaptive  beamforming  stage. 


The  data  parameters  of  the  adaptive  beamforming  stage  are  listed  in  Table  5. 


TABLE  5 

Adaptive  Beamforming  Stage  Parameters 


Name 

Description 

a 

Diagonal  loading  factor 

^bm 

Vector  of  range  gate  indices  for  adaptive  beamforming  stage  training  samples 
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2.4  SUBBAND  PULSE  COMPRESSION 


Input:  N^^i,  subband  data  cubes  of  size  x  x  . 

Output:  subband  data  cubes  of  size  x  x  . 

The  pulse  compression  stage  filters  the  data  to  concentrate  the  signal  energy  of  a  relatively  long 
transmitted  radar  pulse  into  a  relatively  short  pulse  response.  This  compression  is  accomplished  by 
applying  a  FIR  filter  operation  to  the  range  data  for  each  pulse  and  channel  within  each  subband  radar  data 
cube.  That  is,  in  Matlab  notation,  it  operates  on  each  a:  =  C(i,  k,  s),  each  of  which  are  of  length  ;  this 

partitioning  is  illustrated  in  Figure  7.  The  FIR  filter  has  coefficients  I e  {O...A/p^- 1}.  The  output  of 

the  filter,  the  vector  y  of  length  ,  is  the  convolution  of  the  filter  with  the  input  a  : 

yU]  =  X  xU-l]hp^ll],forj  =  0,  1,  ...N^g  .  (6) 

/  =  0 


Subband  s 


Figure  7.  Data  partitioning  of  input  data  cube  for  pulse  compression  stage. 


Unlike  the  time  delay  and  equalization  stage,  the  number  of  taps  in  the  pulse  compression  FIR  filter 
is  usually  much  larger  than  four  or  so  taps.  Therefore,  it  is  usually  more  efficient  to  implement  this 
convolution  with  the  fast  convolution  technique  using  FFTs,  possibly  even  considering  an  overlap-and- 
save  method  [4]. 


13 


The  subband  pulse  compression  stage  data  parameter  is  listed  in  Table  6. 


TABLE  6 


Subband  Pulse  Compression  Stage  Parameters 


Name 

Description 

Number  of  taps  in  the  subband  pulse  compression  filter 

2.5  DOPPLER  FILTERING 

Input:  subband  data  cubes  of  size  N^,„  x  x 

Output:  subband  data  cubes  of  size  x  x  x 

The  Doppler  filter  stage  processes  the  data  so  that  the  radial  velocity  of  targets  relative  to  the 
platform  can  be  determined.  Doppler  filtering  accomplishes  this  using  a  straightforward  FFT  applied  to  the 
pulse  data  for  each  subband  range  gate  and  beam.  A  further  step  made  in  the  Doppler  filtering  which  aids 
the  STAP  stage  in  filtering  out  ground  clutter  involves  taking  staggered  pulse  sample  sets;  is  the 
number  of  pulse  staggers.  Using  g  ( i  <  g  <  as  the  index  of  the  stagger  number,  the  staggers  generated 
in  Matlab  notation  arey^  ^  =  C(i,j,  g'(Np„-N^,^g+g),  s);yij  g  is  a  vector  of  length  For  example,  with 
^siag  =  3  the  first  input  stagger  vector  would  include  the  first  pulse  sample  up  to  the  Afp^,-2nd  pulse 
sample,  the  second  input  stagger  vector  would  include  the  second  pulse  sample  up  to  the  1st  pulse 
sample,  and  the  third  input  stagger  vector  would  include  the  third  pulse  sample  up  to  the  th  (last)  pulse 
sample.  This  is  illustrated  in  Figure  8.  With  A^,^^  =  l,  there  is  only  one  stagger  which  is  the  entire  set  of 
pulse  samples  for  a  given  beam  and  subband  range  gate. 

The  output  of  the  Doppler  filter  is 


for  each  subband,  beam,  subband  range,  and  stagger.  Each  is  a  vector  of  size 
The  parameters  of  the  Doppler  processing  stage  are  listed  in  Table  7. 


(7) 
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Nbm  Beams 


i^gger  1  -samples  i  thru  N^2 

l^lagger  2  ^  sampti^&i  thru 

3  r-  aano^es  3  timi 

r'v 

Subband  s 


Figure  8.  Data  partitioning  showing  three  staggers  of  input  data  cube  for  Doppler  filtering  stage. 


TABLE  7 

Doppler  Processing  Stage  Parameters 


Name 

Description 

^stag 

Number  of  pulse  staggers 

2.6  SPACE-TIME  ADAPTIVE  PROCESSING  (STAP) 

Input:  subband  data  cubes  of  size  x  x  x 

Output:  j  subband  data  cubes  of  size  x  x 

The  STAP  stage  is  a  second  beamforming  stage  which  removes  further  radar  spatially-localized 
interference  and  ground  clutter  (space-  and  time-distributed)  interference.  Using  the  z,  ^  ^  outputs  from  the 
Doppler  filtering  stage, 

'  ■ 

^j.k,g=  -  •  (8) 

The  size  of  Xj  ^  ^  is  x  1,  and  it  is  one  stagger  of  samples.  Next,  the  staggers  are  stacked  upon  each 
other  to  get  a  space-time  snapshot: 
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(9) 


hk  = 


for  the  test  cells,  which  makes  x,.  *  is  We  will  use  these  vectors  of  size  i  as  our 

input  data  to  the  STAP  filter,  but  we  shall  first  build  the  adaptive  filter  matrix.  Defining 


-jlnr 

e  =  0 . AT,;,-! 


(10) 


and  using  the  adaptive  beamforming  matrix  of  the  subband  data  cube  currently  being  processed,  of 
size  N^i^x  define  the  steering  vector  v,  with  length  as: 

»'  =  (<i)\  •  (11) 


The  steering  vector  v  must  be  expanded  to  become  size  x  Refer  to  the  Matlab  code  for 

the  actual  operations;  for  the  sake  of  discussion,  we  can  define  the  x dimensioned  steering 

matrix  V  =  repmat(v,  N^„f).  To  simplify  the  discussion,  let  us  define  =  V(:,  b),  where  1  <  A  <  . 

Now  we  will  choose  the  space-time  snapshots  for  the  training  data  from  this  set: 


Stop 


and  with  it  build  the  interference  estimate  covariance  matrix: 

^k’^  ^  'L\  k^%  k  • 
'(/=  1 

Now  the  corresponding  adaptive  weights  are  calculated  as: 


(12) 


(13) 


'k,b 


=  r: 


-1 

t  ^b 


(14) 


As  in  adaptive  beamforming,  the  above  equation  is  referred  to  as  the  Wiener-Hopf  equation  for 
calculating  adaptive  filter  weights.  While  this  explains  the  mathematics  behind  calculating  the 
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interference-nulling  matrix,  calculating  the  inverse  of  is  not  used.  Instead  LQ-decomposition  is  used  to 
calculate  each  of  the  ^  vectors. 

As  explained  in  Appendix  A,  calculate  =  IqhouseiR/^);  is  a  lower  triangular  f^bm^stag^^bm^stag 
matrix.  Now  using  t  =  ifw*  *  as  a  temporary  vector  of  size  x  l ,  we  can  solve  L^i  =  vj  for  t  using 

forward-substitution  and  =  t  for  ^  using  back-substitution.  Solving  this  for  every  b  =  \...Ncnb 

provides  ^  ...  for  Doppler  bank  k',  is  a  matrix  of  size  x  Of  course  all  of 

these  calculations  must  be  done  for  each  Doppler  bank,  k  = 

Refer  to  Appendix  A  for  an  explanation  of  this  method  of  solving  the  Wiener-Hopf  equation. 

Finally,  the  adaptive  weights  vectors  are  applied  to  the  snapshots  jc^.  ^  of  size  x  l; 


yj,k 


—  W^Xji^ 


(15) 


The  resulting  vector  yj  ^  is  of  size  N^^b  ^  l  must  be  calculated  for  each  J=l-Nsrg  and 

-f^dop- 

The  STAP  stage  data  parameters  are  listed  in  Table  8. 

TABLE  8 


Space-Time  Adaptive  Processing  Stage  Parameters 


Name 

Description 

M 

^  ''stag 

Number  of  pulse  staggers 

Ks 

Number  of  training  data  set  vectors 

T 

Stag 

Vector  of  space-time  snapshot  indices  for  the  training  data  set  for  space-time 
adaptive  processing 

2.7  SUBBAND  SYNTHESIS  (COMBINING) 

Input:  Narrowband  subband  data  cubes  of  size  x  Af^^^x 

Output:  Processed  wideband  radar  data  cube  of  size  N^„b  ^  Kg  ^  Kop- 

Since  the  subband  analysis  and  synthesis  are  signal  processing  duals  of  one  another,  the  subband 
synthesis  stage  was  covered  in  Section  2.1. 
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2.8  TARGET  DETECTION 


Input:  Processed  radar  data  cube  of  size  N^„i,  x  x 

Output:  Coordinates  of  detected  targets  in  radar  data  cube. 

Though  there  are  many  techniques  for  detecting  targets  in  the  processed  radar  data  cube  including 
various  integrators  and  hard  limiters,  this  application  shall  use  a  version  of  constant  false  alarm  rate 
(CFAR)  detection  with  three-dimensional  grouping.  The  data  parameters  of  the  target  detection  stage  are 
listed  in  Table  9. 

2.8.1  CFAR  Detection 

During  CFAR  detection,  a  local  noise  estimate  is  computed  from  the  range  gates  near  the  cell 
C(i.j,  k)  under  test.  A  number  of  guard  gates  G  immediately  next  to  the  cell  under  test  will  not  be  included 
in  the  local  noise  estimate  (this  number  does  not  affect  the  throughput).  The  range  cells  involved  in 
calculating  the  noise  estimate  for  a  particular  Doppler  bin  and  beam  are  shown  in  Figure  9.  For  each  cell 
C(i,j,  k),  the  value  of  the  noise  estimate  T(i,j,  k)  is  calculated  as: 

ni,ik)  =  ,^  y  \cii,j^i,kf+\c(i,j-i,k)\^  .  (16) 

^  ^cfar 

/  =  C+1 


Beams 


C(i,i.k) 

Nc„  G  I  G  N. 


T(i,j.k) 


t 


Figure  9.  Data  partitioning  of  input  for  CFAR  detection  stage. 
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At  the  boundaries  when  j-{G  +  - 1 )  <  1  or  y  +  ( G  +  >  N,g,  only  the  side  range  gate  cells  that 

are  within  the  data  cube  are  included,  and  the  denominator  of  the  leading  scaling  factor  is  decreased  to 
reflect  the  number  of  side  range  gate  cells  that  were  actually  used  for  the  calculation.  For  each  cell  C(i,  j, 
k),  the  quantity  \C{i,j,  k)\^/T{i,J,  k)  is  calculated;  this  represents  the  normalized  power  in  the  cell  under  test. 
If  this  normalized  power  exceeds  a  threshold  p  ,  the  cell  is  considered  to  contain  a  target.  This  threshold 
can  be  adaptively  determined,  but  for  this  application,  we  shall  set  it  to  a  constant.  The  only  cells  in  the 
data  cube  that  are  processed  by  the  following  stages  are  those  where  a  target  has  been  detected. 

2.8.2  Three-dimensional  Grouping 

During  this  step,  detection  in  adjacent  cells  are  grouped  to  avoid  having  multiple  detections 
associated  with  the  same  target.  The  power  of  each  target  detected  by  the  CFAR  algorithm  will  be 
compared  to  the  power  of  each  cell  in  the  3x3x3  cube  centered  on  that  cell.  Therefore,  each  cell  under 
test  will  be  compared  to  26  other  cells  (27  -  1 ,  where  27  =  3^  and  1  cell  is  the  cell  under  test,  the  cell  in  the 
middle  of  the  cube).  If  the  detection’s  raw  power  is  the  local  maximum,  it  is  retained.  Otherwise,  the 
detection  is  discarded;  it  is  considered  grouped  with  the  detection  in  the  cell  with  the  higher  power. 

TABLE  9 


Target  Detection  Stage  Parameters 


Name 

Description 

^cfar 

Number  of  range  gates  used  in  forming  the  noise  estimate  from  each  side  of 
the  cell  under  test 

G 

Number  of  guard  cells  on  both  sides  of  the  cell  under  test 

Normalized  target  power  threshold 

2.9  TARGET  PARAMETER  ESTIMATION 

Input:  Coordinates  of  detected  targets  in  radar  data  cube. 

Output:  A  target  report. 

This  stage  of  the  GMTI  processing  packages  the  target  parameters  for  the  target  report.  As  presented 
in  Section  1.4,  a  target  report  consists  of  the  following  data  for  each  target: 

•  target  power  or  signal-to-noise  ratio  (SNR) 

•  target  azimuth  (in  degrees  or  radians) 

•  target  range  (in  meters  or  kilometers) 

•  target  radial  velocity  (in  m/s  or  km/h) 
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♦  target  time  stamp 

The  target  power  is  the  quantity  \C(i,j,  k)^ k)  for  each  target  cell  C(i,j,  k)  which  was  calculated 
by  the  CFAR  detection.  The  target  azimuth  is  determined  by  the  beam  index  number  and  the  platform 
position  and  heading.  We  shall  assume  a  temporarily  stationary  platform,  so  it  will  suffice  to  just  use  the 
beam  index  number.  The  target  range  is  determined  by  the  range  gate  index  number  and  the  platform 
altitude,  grazing  angle,  and  azimuth  angle.  Again,  we  shall  assume  a  temporarily  stationary  platform,  so  it 
will  suffice  to  just  use  the  range  gate  index  number.  The  target  radial  velocity  is  determined  by  the  Doppler 
index  number  and  the  speed  and  heading  of  the  platform,  but  with  the  assumed  temporarily  stationary 
platform,  passing  the  Doppler  index  number  will  suffice.  The  target  time  stamp  will  be  reported  as  the 
index  to  the  data  cube  from  which  this  target  was  determined. 
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3.  MATLAB  IMPLEMENTATION 


The  GMTI  processing  chain  will  be  initialized  with  a  fiinction  with  the  following  signature: 
GMTI_params  =  GMTI_initialize ( ) ; 

This  function  will  load  control  parameters  outlined  in  Section  1 .2  into  the  structure  GMTI_params. 
The  GMTI  processing  chain  will  be  invoked  with  a  function  with  the  following  signature: 

targets  =  GMTI (input_data,  GMTI_params) ; 

This  function  call  passes  the  aforementioned  GMTI_params  and  the  input_data  outlined  in 
Section  1 .3,  and  it  returns  the  target  structure  outlined  in  Section  1 .4. 

In  order  to  keep  memory  usage  to  a  reasonable  level,  the  GMTI  function  will  process  only  one  radar 
data  cube  at  a  time,  and  within  that  radar  data  cube,  it  will  process  only  one  subband  radar  data  cube  at  a 
time.  Each  of  the  stages  depicted  in  Figure  1  will  have  its  own  function  call  within  the  GMTI  function. 


APPENDEXA 

SOLVING  THE  WIENER-HOPF  EQUATION 


The  Wiener-Hopf  equation  is  used  in  both  the  adaptive  beamforming  and  STAP  stages.  The  equation 
w  =  defines  an  adaptive  filter  vector,  w,  as  the  product  of  the  inverted  estimated  noise  covariance 
matrix,  R ,  and  the  nonadaptive  steering  vector,  v .  ^  is  composed  using  a  set  of  training  vectors,  xj.,  such 
that  R  =  . 

However,  R  and  its  inverse  are  seldomly  used  to  solve  the  Wiener-Hopf  equation  [3]  because 
calculating  R  and  its  inverse  causes  numerical  instability.  Instead,  LQ-decomposition  with  forward-  and 
back-substitution  is  used.  First,  the  more  general  solution  of  the  STAP  Wiener-Hopf  equation  is  discussed, 
and  then  the  adaptive  beamforming  solution  (in  which  a  diagonal  loading  factor  is  added)  is  shown. 

Compute  the  LQ  decomposition  xj.  =  ig ,  where  L  becomes  a  lower  triangular  matrix  and  Q 
becomes  a  unitary  matrix.  The  LQ  decomposition  is  related  to  the  QR  decomposition  described  in  [2],  by 
inputting  the  Hermitian  of  Xj.,  and  receiving  the  output  as  the  Hermitian  of  L.  Matlab  code  for  the  LQ 
decomposition  follows: 

%lqHouse  Perform  LQ  decomposition  of  input  matrix 
function  L  =  IqHouse (A) 

%  Find  the  size  of  the  input  matrix  A. 

[m,  n]  =  size  (A)  ; 

%  Loop  over  each  row  of  A. 
for  row  =  l:m 


%  Compute  the  Householder  vector  v. 

X  =  0  ; 

x(row:n,  1)  =  (A(row,  row:n) ) ' ; 

V  =  x; 

v(row)  =  v(row)  +  (x(row)  /  abs(x(row)))  *  norm(x) ; 
beta  =  -  2  /  (v'  *  v)  ; 


%  Apply  the 
matrix. 

w ( row : m , 
A (row :m, 


Householder  reflection  to  the  remainder  of  the 

1)  =  beta  *  A(row:m,  row:n)  *  v(row:n) ; 

row:n)  =  A(row:m,  row:n)  +  (w(row:m)  *  v{row:n)'); 
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end  %  for  row 


L  =  A; 

Now  the  estimated  noise  covariance  matrix  can  be  expressed  as: 


R  =  =  (LQ){LQf  =  LQ(Q^l'^)  =  LZ,"  ,  (17) 


and  the  Wiener-Hopf  equation  can  be  rewritten  as: 

w  =  R  V  =  {LL^)  *v 
(,LL^)w  =  =  V 


(18) 

(19) 


Letting  t  =  L^w,  Lt  =  v  can  be  solved  for  t  by  using  forward  substitution  since  L  is  lower  triangular. 
Then  using  t  =  l'^w  ,  a  solution  for  the  adaptive  filter  vector,  w ,  can  be  found  by  using  back-substitution 
since  is  upper  triangular. 

Now  the  added  complication  of  the  added  diagonal  loading  factor  in  the  adaptive  beamforming  stage 
is  discussed^To  solve  this  problem,  define:  x/  =  .  Now  x-,!  =  L’Q ,  and  the  Wiener-Hopf  equation, 

M>  =  ((L')(L')  )  V ,  can  be  solved  as  described  above. 

Finally,  the  denominator  of  equation  (4)  is  a  scalar  which  can  be  solved  in  the  following  manner. 
Using  R  =  (L'XL'f, 


v^r'v  =  v%L')iL’fi'v  = 


(20) 


Letting  u  =  (U)  v,  L'u  =  v,  which  can  be  solved  for  w  using  forward  substitution,  and 

H  c  7 

V  R  V  =  u  u,  which  is  a  scalar. 
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ACRONYMS 


CFAR  - 
CPI  - 
FFT  - 
FIR  - 
GMTI  - 
PCA  - 
PRI  - 
SNR  - 
STAP  - 


Constant  False  Alarm  Rate 
Coherent  Processing  Interval 
Fast  Fourier  Transform 
Finite  Impulse  Response 
Ground  Moving  Target  Indicator 
Polymorphous  Computer  Architecture 
Pulse  Repetition  Interval 
Signal-to-Noise  Ratio 
Space-Time  Adaptive  Processing 
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